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Resumen 


Se muestra una aplicación de los polinomios splines cúbicos ala resolución numérica 
de la ecuación diferencial 1D de la dispersión. Para el caso no estacionario, la 
variación en el tiempo se refleja en los coeficientes del polinomio spline y la 
discretización en el tiempo se hace aplicando el Método de las Diferencias Finitas. 
Todos los algoritmos son implementados computacionalmente en el asistente 
matemático Matlab. Se comparan los resultados obtenidos de las soluciones 
aproximadas con las mediciones de un caso real de una prueba de trazador y 
además con la solución del Método de los Elementos Finitos. 
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Introducción 


Los polinomios splines tienen gran aplicación en 
la interpolación de funciones. En el problema 
de la interpolación con los polinomios splines es 
necesario añadir condiciones adicionales en los 
extremos del intervalo. Samarski (1986) plantea 
el problema de la interpolación por tramos 
1D a partir de los polinomios splines cúbicos 
de forma que debe resolver un sistema lineal 
de ecuaciones con matriz tridiagonal. Para 
obtener solución única añade dos condiciones 
adicionales, asociadas con la segunda derivada 
en los extremos, lo cual se corresponde con 
el llamado spline natural (Samarski, 1986). 
Álvarez et al. (2004) plantean la obtención 
del polinomio cúbico, integrando dos veces 
un polinomio lineal que interpola los valores 
de la segunda derivada en un intervalo. Los 
splines también se aplican en la resolución de 
ecuaciones diferenciales, con lo cual se obtiene 
que las funciones solución cumplen condiciones 
de continuidad de las derivadas en puntos 
interiores de la región sin necesidad de tener 
que conocer de forma previa los valores de 


las derivadas en dichos puntos. Szymkiewicz 
(1993) descompone la ecuación 1D del trans- 
porte en dos ecuaciones, una con la parte 
advectiva y otra ecuación con la parte difusiva 
(Szymkiewicz, 1993). Utiliza el método de las 
características junto con los splines cúbicos para 
resolver la ecuación de la parte advectiva y 
utiliza el Método de los Elementos Finitos para 
resolver la ecuación de la parte difusiva. Tsai 
et al. (2004) también utilizan los polinomios 
splines cúbicos para resolver la ecuación 1D 
no estacionaria del transporte (Tsai et al., 2004) 
y emplean diferentes tipos de condiciones 
adicionales en los extremos del intervalo. 

El objetivo de este trabajo es mostrar la 
resolución de la ecuación diferencial 1D de la 
dispersión, utilizando los splines cúbicos para un 
caso real de estudio de una prueba de trazador, 
así como evidenciar las ventajas de aplicar este 
método que fundamentalmente radica en que la 
solución obtenida por este método garantiza la 
suavidad de la misma en cada uno de los puntos 
de la discretización de la región de estudio, lo 
cual no ocurre con el Método de los Elementos 
Finitos. 
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La función spline cúbica 


Se tiene una red formada por los subintervalos 
[x, x,,], con i = 1,...N,, de un intervalo [0, L] 
La función u = u(x) definida en dicho intervalo 
se puede aproximar por un spline cúbico, 
expresado en función del extremo izquierdo 
del intervalo de la forma: 


9/0) =0 +0 (x-x,)+ a) (x-x,) 
+a(x-x,) para ¡=1,...,N, (1) 


Las condiciones de interpolación que debe 
cumplir el spline son las siguientes: 


p(x,)=4, parai=1,..,N, (2) 
Pm, (Xn 11)=n +1 (3) 


Las condiciones de continuidad del spline 


cúbico y las condiciones de suavidad de su 
primera y segunda derivadas son las siguientes: 


di(X;,1) = dia %;,1) parai=1,...,N,-1 (4) 
9'i(%;,1) =isi[X;,1) parai=1,...,N,-1 (5) 
0% (1) = Oi (%;,1) parai=1,...,N,-1 (6) 


Como puede observarse de (2) y de (3) se 
tiene N, + 1 condiciones de interpolación, y 
de (4) hasta (6) se tiene 3(N, - 1) condiciones 
de continuidad; por lo tanto se tienen 4N_ - 2 
condiciones en total y hay 4N_ coeficientes 
incógnitas. Teniendo en cuenta las condiciones 
de (1) a (6) se obtiene el siguiente sistema de 
ecuaciones: 


i 2 3, 

Ax a,+Ax” a3+Ax" a, =84; 4; 
i i 2.5 _ ir 
4,+2 Ax a3+3 Ax” a, =4) 

i i ¡+1 
43+3 Ax a, =43 


NN 
Axa, +Ax* ay +Ax ay =Uy Uy 0) 
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donde¡=1,..., N_—1 y Ax representa la longitud 
del intervalo. 

En el problema de la interpolación se 
conocen todos los valores u, con ¡=1,..., N + 
1, de la función a interpolar u = u(x) en los N, 
+ 1 nodos. Los valores de los coeficientes a; 
correspondientes a los términos independientes 
de (1), coinciden con los valores de la función 
incógnita u; por lo tanto, el número de 
coeficientes incógnitas se reduce a 3N_ y el 
número total de ecuaciones se reduce a 3N, - 2. 
Se le añaden dos condiciones adicionales en los 
extremos para que el problema tenga solución 
1 1 
las tres primeras ecuaciones del sistema (7) 


única. Para cada valor de í, con i¡=1...., N 


Ed 


se pueden escribir en notación matricial en la 
forma: 


a 


donde C y B son matrices, y A y E son vectores, 
según se muestran a continuación: 


Ax Ax? Ax? 00.0 
C=| 1 2Ax 3Ax? B=|1.0:0 
0 1  3Ax 01.0 
a 44174 
A=4 a F=i 0 
d 0 


La resolución del sistema de ecuaciones (8), 
en conjunto con las dos condiciones adiciona- 
les en los extremos del intervalo, da la solución 
del problema de interpolación. 


Aplicación de los polinomios splines 
cúbicos a la resolución de la ecuación de 
la dispersión 1D no estacionaria en un 
medio poroso 


Se considera la ecuación de la dispersión 1D 
no estacionaria en un medio según la forma 
(Marón, 2002): 


Marón-Domínguez y Gutiérrez-de-la-Rosa, Resolución de la ecuación de la advección-dispersión en una dimensión, usando... 


2 
a y EA A 
Ox ax? ot 


para0<x<L, t>0 


(9) 
donde O es la porosidad del medio; R,, el 
coeficiente de retardo; C = C (x, £), la función 
incógnita que representa la concentración; 
U, la velocidad del agua; D, el coeficiente de 
dispersión, y A es el coeficiente de degradación 
o de decaimiento. A la ecuación anterior 
hay que añadirle una condición inicial y dos 
condiciones de frontera: 
C(0,t)=Ce(t)  C(L,t)=Cs(t)  C(x,0)=CUx) 
Como este problema es no estacionario, 
entonces la variación en el tiempo se propone 
que sea recogida en los coeficientes del spline 
cúbico, de forma que el polinomio tiene la 
forma: 


9 (1,1) =0,(0)+03(t) (x-x,)+a3(+) lua]? 


+0,(t) Ez (10) 
Sustituyendo (10) en (9), para X =X, se 
obtiene la expresión: 


0Rr (1+ At Aa" +(At 0 u)a,' +(-2 At D)az' 


=0 Rya¡? (11) 


donde ¿ =2,..., N, y k=2,..., N, + 1. Para 
el polinomio del primer tramo del spline 
correspondiente a ¡= 1 se obtiene: 


(At 15) U) a, +(-2 At D) a 


=0Rpa "0 Rp(1+AtA) CI (12) 
la expresión anterior tiene incluida la condición 
de frontera en el extremo izquierdo. Evaluando 
el polinomio spline del último tramo en X = 


Xw,17 Se obtiene: 


OR, (1+AtA) a F+(At8 u) ads k 

+ (22 At D) ab: “0 Rae pel 

Aaa A? ad acia 
OR (1+2 At) ae" 

+(0 R¿Ax+ AU +AtAO RpAx) a) *" 
+(2 Ax AtO0 U-2 44D +0 Ryax? 


+ ARO Rad Jaj" d (13) 


+(3 AY'AtO U-6 Ar A£D+0 Rar? 
+Ath 9 RyAx*) ayek 


N, k-1 N, kl N, k-1 N, k1 
= ORp (aj +Axa),* + Ax? a: +A as” ] 


Las ecuaciones (11), (12) y (13) expresan 
que se debe satisfacer la ecuación diferencial 
(9) en todos los puntos de la red. Para cada 
tiempo, el sistema final de ecuaciones se 
obtiene ensamblando el sistema (8), junto con 
las ecuaciones (11), (12) y (13). 

Para el cálculo de la solución en el primer 
intervalo de tiempo es necesario calcular los 
valores de los coeficientes de los splines en 
todos los tramos. Para esto se hace uso de la 
condición inicial, utilizando la interpolación 
propuesta por el sistema (8). 


Comprobación de los algoritmos 
propuestos utilizando una aplicación 
práctica 


Con el objetivo de comprobar los algoritmos 
propuestos en este trabajo se hace uso de 
los mismos datos o mediciones reales de la 
prueba de trazador detallada en el artículo 
de Marón (2002). En el citado trabajo, los 
datos fueron utilizados para la calibración 
o estimación de los parámetros hidráulicos 
de un humedal construido en el laboratorio. 
Marón utilizó, para la discretización en 
el espacio, el Método de los Elementos 
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Finitos (MEF), con ayuda de los elementos 
cuadráticos y cúbicos. Esto permite que la 
comparación del método desarrollado en este 
trabajo con los splines se haga no sólo con los 
datos reales observados sino también con 
los resultados obtenidos por el MEF cúbico, 
ya que los polinomios utilizados para la 
construcción del spline fueron cúbicos. Para 
la discretización en el tiempo se hizo uso del 
Método de las Diferencias Finitas, al igual que 
en el artículo citado. El algoritmo de cálculo 
propuesto fue implementado en el asistente 
matemático Matlab, el cual se acopló con un 
algoritmo de optimización, que ayudó a la 
estimación de los parámetros de dispersión, 
retardo y degradación (D, R, y A). 

En el cuadro 1 se muestran los resultados 
de la calibración de dichos parámetros con 
el algoritmo de los splines y se comparan con 
los resultados obtenidos por el MEF cúbico, 
obtenidos por Marón en su trabajo (Marón, 
2002). 

En el cuadro 1 se aprecia la no existencia 
de diferencias sustanciales entre los valores 
obtenidos por el MEF desarrollado por Marón 
(Marón, 2002) y el método utilizado por los 
autores. Se observa que no se detectó retardo en 
ambos casos, y se aprecia buena aproximación 
en los resultados de las estimaciones de los 
coeficientes de dispersión y degradación. 

En la figura 1 se muestran las concentra- 
ciones observadas en las secciones x = 0.9 m 
y x= 1.9 m, correspondientes a los puntos de 
observación, en conjunto con las obtenidas 
según el algoritmo desarrollado en este trabajo 
a partir de los splines con polinomios cúbicos 
y las obtenidas según el MEF con elementos 
cúbicos, referidas por Marón en su trabajo. 

En la figura 1 se observa que se han 
obtenido buenas aproximaciones entre las 
concentraciones observadas, (Cobs(*)), y las 
(Ccal(o)), 
(gráficas de la derecha). Es bueno señalar, 


calculadas, en ambas secciones 
además, que la comparación se puede hacer 
entre los valores calculados con los splines y 
los valores obtenidos por el MEF referidos 


Cuadro 1. Valores de los parámetros calibrados. 


MEF (cúbicos) Splines (cúbicos) 
D (m?/min) 0.00025588 0.00018295 
R, 1 1 
2 (1/min) 0.0003008 0.000422 


en el trabajo de Marón (Marón, 2002). Para el 
algoritmo correspondiente a los splines cúbicos, 
se calculó el número de Peclet, cuyo valor fue 
2.5571. 


Conclusiones 


Para el ejemplo real utilizado se evidenció que 
se obtuvieron resultados satisfactorios para 
valores del número de Peclet mayor que 1. 

El método desarrollado por los splines 
cúbicos aporta resultados superiores al Método 
de los Elementos Finitos, ya que la función 
solución satisface en los nodos interiores de la 
red condiciones de suavidad hasta la segunda 
derivada, lo cual no se cumple en el Método de 
los Elementos Finitos. 

Queda esbozado el uso de un método que 
muestra resultados satisfactorios para que un 
especialista en el estudio de aguas residuales 
disponga de una herramienta más. 
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Figura 1. Gráficas de las concentraciones calculadas y las observadas en x 
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Abstract 


MARÓN-DOMÍNGUEZ, D.E. € GUTIÉRREZ-DE-LA-ROSA, A. Solving an one- 
dimensional advection-dispersion equation with spline functions. Water Technology and 
Sciences (in Spanish). Vol. IV, No. 4, September-October, 2013, pp. 173-178. 


The application of cubic spline polynomials to the numeric solution of a one-dimensional 
differential dispersion equation is presented. For the non-stationary case, the variation in time 
was reflected by the spline polynomial coefficients and the discretization in time was obtained 
using the finite difference method. All the algorithms were implemented with the MATLAB 
mathematical assistant. The results were compared to tracer test measurements from a real 
case as well as with the finite element method solution. 


Keywords: cubic splines, interpolation, differential equation, dispersion. 
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